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Abstract —Tht mass redistribution in the earth as a result of an earthquake faulting changes the 
earth's inertia tensor, and hence its rotation. Using the complete formulae developed by Chao and 
Gross (1987) based on the normal mode theory, we calculated the earthquake-induced polar motion 
e.xcitation for the largest 1 1,015 earthquakes that occurred during 1977,0 1993.6. The seismic excitations 
in this period are found to be two orders of magnitude below the detection threshold even with today s 
high precision earth rotation measurements. However, it was calculated that an earthquake of only one 
tenth the size of the great 1960 Chile event, if happened today, could be comfortably detected in polar 
motion observations. Furthermore, collectively these seismic excitations have a strong statistical ten- 
dency to nudge the pole tow'ards 140 F., away from the actually observed polar drift direction. This 
non-random behavior, similarly found in other earthquake-induced changes in earth rotation and 
low-degree gravitational field by Chao and Gross ( 1987), manifests some geodynamic behavior yet to 
be explored. 
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/. Introduction 

The earth’s rotation varies slightly with time. The 3-D earth rotation variation 
can be conveniently separated into two components: (i) The 1-D variation in the 
spin rate, often expressed in terms of the length-of-day variation, (ii) 2-D variation 
in the rotational axis orientation, generically called the nutation when viewed from 
the inertial reference frame, or the polar motion as seen in the terrestrial reference 
frame. 

There are two dynamic types of earth rotation variations (Munk and Mac- 
Donald, 1960). (i) The “astronomical variations” due to external luni-solar tidal 
torques that change the earth’s angular momentum: Well-known examples include 
the tidal braking that causes the earth’s spin to slow down over geological times, 
and the astronomical precession and nutation caused by torques exerted on the 
earth’s equatorial bulge, (ii) The “geophysical variations” caused by large-scale 
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mass movement of internal geophysical processes under the conservation of angular 
momentum. These include tidal deformations in the solid earth and oceans, 
atmospheric fluctuations, hydrological variations, ocean currents, earthquake dislo- 
cations, post-glacial rebound, mantle tectonic movement, and core activities. 

The present paper deals with the polar motion and its seismic excitation. 
Possible interactions between seismicity and earth rotation have been under consid- 
eration since the discovery of the polar motion a hundred years ago (e.g., 
Lambeck, 1980). On one hand, a non-uniform earth rotation would give rise to a 
time-varying stress field inside the earth, which in turn may affect the triggering 
process of earthquakes. This effect with respect to the polar motion is in general an 
order of magnitude smaller than its tidal counterpart (Lambert, 1925), but having 
a longer time scale on the order of a year. Chao and Iz (1992) conducted a 
statistical test in an attempt to correlate the occurrence of earthquakes with the 
contemporary polar motion based on nearly 10,000 earthquakes that occurred 
during 1977-1991; but only a weak correlation was found. A more definitive 
conclusion awaits further studies taking into account the tensorial nature of the 
seismic source mechanism and the stress field induced by polar motion. 

On the other hand, the mass redistribution in the earth as a result of an 
earthquake faulting changes the earth’s inertia tensor, and hence its rotation. Early 
simplistic earthquake faulting models considering only regional dislocations greatly 
underestimated the effect (Munk and MacDonald, 1960). After the 1964 Alaska 
earthquake it was recognized that the seismic dislocation, although decreasing with 
focal distance, remained non-zero even at teleseismic distances away from the fault 
(Press, 1965). When integrated globally (see equation (2) below) considering the 
size and source mechanism of the earthquake this displacement field can give finite 
effects (Mansinha and Smylie, 1967). Smylie and Mansinha (1971), Dahlen 
(1971, 1973), Rice and Chinnery (1972), and O’Conneli. and Dziewonski 
(1976), among others, calculated the seismic excitation of polar motion based on 
realistic models. It was concluded that the great i960 Chile and 1964 Alaska 
earthquakes should have produced significant changes in polar motion (see below). 
Unfortunately the accuracy of the polar motion record at that time was insufficient 
to yield any conclusive detection. 

In the last three decades, the measurement of earth rotation has been progres- 
sively improved by three orders of magnitude in both accuracy and temporal 
resolution (e.g., Eubanks, 1993, for a review). This is achieved through the 
advances in modern space geodetic techniques, primarily Very-Long-Baseline Inter- 
ferometry (VLBI), Satellite Laser Ranging (SLR), and Global Positioning System 
(GPS). In fact, these techniques have measured the contemporary relative tectonic 
plate motions to within a few mm/year, and revealed coseismic, near-field displace- 
ments caused by recent earthquakes. The SLR technique has also greatly benefitted 
the determination of the earth’s global gravitational field and detected slight tem- 
poral variations in the low-degree field. In the area of earth rotation measure- 
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ment, the current accuracy is estimated to be within 1 milliarcsecond ( 1 mas = 
4.85 X 10“"^ radian, corresponding to about 3 cm of distance on the surface of the 
earth), while the formal errors are as low as 200 /t as. 

Taking advantage of these modern data as well as the seismic centroid moment 
tensor solutions made available in the Harvard CMT catalog (for a review, see 
DziEWONSKi et ai. 1993), Souriau and Cazenave (1985) and Gross (1986) 
calculated the seismic excitation of polar motion using Dahi.en'S (1973) formula- 
tion for post- 1977 major earthquakes. It was concluded that the polar motion 
excited by those earthquakes were all too small to be detected even in the modern 
earth rotation data. Chao and Gross (1987) developed complete formulae for 
calculating the earthquake-induced changes in earth's rotation and low-degree 
gravitational field based on the normal mode theory of Gilbert ( 1970). They made 
calculations for the 2,146 major earthquakes that occurred during 1977-1985. In 
addition to confirming the above conclusion, the result indicated that earthquakes 
have a strong statistical tendency to make the earth rounder and more compact. It 
also demonstrated a similar tendency for earthquakes to nudge the rotation pole 
towards the direction of about 150 E. Chao and Gross (1995) and Chao et aL 
( 1995) extended the calculation to include 1 1,015 major earthquakes for the period 
1977-1993 to evaluate the corresponding rotational and gravitational energy 
changes. 

The present paper will study the seismic excitation of polar motion in the 
framework of Chao and Gross (1987) but making use of the updated result of 
Chao and Gross ( 1995). The corresponding changes in length-of-day will only be 
presented in passing, as they are in general two orders of magnitude smaller (that 
is to say, polar motion is hundreds of times easier to excite than lenglh-of-day by 
earthquakes). The physical reason is the following. The geophysical excitation acts 
against the “inertia" of the earth system. For length-of-day the inertia is the earth's 
axial moment of inertia C, whereas the inertia for polar motion is the difference 
between the axial and equatorial moments of inertia, C -A (see equation 2), which 
is only 1/300 of C. 


2. Dynamics ami Calculation 


The excitation of the polar motion is governed by the conservation of angular 
momentum. The equation of motion is customarily expressed as (Munk and 
MacDonald, 1960): 


rn + 


/ dm 
(7 dt 




( 1 ) 


where m is the complex-valued pole location measured in radian; its real part is the 
component (along the Greenwich Meridian) and the imaginary part the v' 
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Polar motion & Excitation 



Year 
Figure I 

The X and y components of the observed polar motion and its excitation function (V', obtained by 
deconvolution) for 1976.4 1994.0, in units of milliarcscconds. The straight lines are the least-squares fit 

to the excitation function. 


component (along the 90 E longitude), a is the frequency of the free Chandler 
wobble with a nominal period of 435 days and a Q value of 100, and 'E is the 
complex-valued excitation function. Mechanically equation (1) is analogous to the 
excitation of a simple harmonic oscillator with a natural frequency a. 

The polar motion m traces out a prograde, quasi-circular path on the order of 
10 m in the vicinity of the North Pole. Figure 1 shows the .v and v positions of the 
pole at nominal daily intervals according to the ''Space93” dataset (Gross, 1994) 
derived from space geodetic measurements during 1976.4-1994.0. Besides a slow 
polar drift, the oscillation consists mainly of the annual wobble and the Chandler 
wobble. It is continually excited (otherwise it would decay away in a matter of 
decades); and the excitation function T can be obtained numerically according to 
equation (1) in a process of deconvolution. The T' thus obtained is also given in 
Figure 1, together with the least-squares fitted straight lines representing the polar 
drift (see below). The geophysical problem is to identify and understand the sources 
of this “observed"’ T. It is now known that a major source is the variation of the 
atmospheric angular momentum (e.g., Chao, 1993; Kuehne et al., 1993). The 
problem is far from closed, and the earthquake dislocation remains a candidate 
excitation source. 
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The polar-motion excitation 4^ due to mass redistribution is given by 

4^= 1.61(A/,, +/ A/,,)/(C-^) (2) 

where / denotes the inertia tensor, the factor 1.61 takes into account the earth's 
non-rigidity and the decoupling of the fluid core from the mantle in the excitation 
process. Note that T should also include an additional term due to mass motion, 
but that term is negligible in the case of abrupt seismic sources (Chao, 1984). 

With an abrupt step-function time history (compared to the considerably longer 
time scale of the polar motion), an earthquake faulting generates a co-seismic, 
step-function displacement fleld u (after the seismic waves have died away). 
Knowing the seismic moment tensor, u anywhere in the earth can be evaluated by 
the normal mode summation scheme of Gilbert (1970). The task, then, is to 
calculate the seismic T according to equation (2), which consists of evaluating A/ 
by a properly weighted integration of u over the globe. The reader is referred to 
Chao and Gross ( 1987) for details of the formulation and calculation method. 
This normal mode scheme is found to be extremely eflflcient. 

Calculation has been conducted for 11,015 major earthquakes (with nominal 
magnitude greater than 5.0) that occurred during 1977.0-1993.6, using the seismic 



Year 


Figure 2 

The A and i components of the seismic excitation of polar motion. The straight lines are the 
least-squares lit to the c.xcilation function. 
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centroid moment tensor solutions published in the Harvard CMT catalog (e.g., 
Dziewonski et aL, 1993). Smaller earthquakes, although numerous, have com- 
pletely negligible effects. The adopted normal mode eigenfrequencies and eigenfunc- 
tions belong to the spherically symmetric earth model 1066B of Gilbert and 
Dziewonski (1975). The net effect is then the accumulation of individual step- 
function contributions: 4^(0 = — /„), (« = 1, . . . , 1 1,015). Chao and 

Gross ( 1995) have also calculated the corresponding changes in length-of-day and 
low-degree gravitational field. 


3. Results and Analysis 


Figure 2 shows the calculated polar-motion excitation function ^(t) by the 
earthquakes. The starting value is arbitrarily chosen to be zero. Comparing Figure 
2 with 1, it is obvious that the magnitude of the seismic excitation is insignificant 
during 1977.0-1993.6: It is of two orders of magnitude too small to explain the 
observed polar-motion excitation. 

For the purpose of illustration. Table 1 lists the results for the following eight 
largest earthquakes in recent decades (with seismic moment A/o exceeding 10^' Nm). 
These results have been reported elsewhere by Chao and Gross (1995). 


Event I: May 22, 1960, Chile 

Event II: March 28, 1964, Alaska, U.S.A. 

Event III: August 19, 1977, Sumba, Indonesia 

Event IV: March 3, 1985, Chile 

Event V: September 19, 1985, Mexico 

Event VI: May 23, 1989, Macquarie Ridge 

Event VII: June 9, 1994, Bolivia 

Event VIII: October 4, 1994, Kuril Is., Russia 


The source mechanism of Events I and II, which occurred before the span of the 
Harvard catalog, are taken from Kanamori and Cipar (1974) and Kanamori 


Table 1 


Magnitude and direction of polar -motion excitation by eight great earthquakes, and the corresponding 

length-of-day changes 


Event 

(A/„, 10-* Nm) 

I 

(270) 

II 

(75) 

HI 

(3.6) 

IV 

(10) 

v 

(1.1) 

VI 

(1-4) 

vn 

(2.6) 

VIII 

(3.9) 

ALOD (/4s) 

-8.4 

6.8 

0.33 

-0.10 

-0.089 

-0.059 

0.192 

-0.053 

I't'l (mas) 

22.6 

7.5 

0.21 

0.18 

0.084 

0.114 

0.331 

0.256 

arg(^) C E) 

115 

198 

160 

no 

277 

323 

122 

129 
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(1970), respectively. Events VII (a deep-focused event ) and VIII in 1994 are also 
outside our studied period. They have the largest seismic moment since Event HI in 
1977. 

The following facts are observed: (i) The earthquake-induced ALOD is small: 
even the largest earthquakes (Events I and II) would only cause a ALOD that is 
barely discernable even with today’s measurement precision, (ii) The direction of 
polar motion excitation, arg(T), tend to cluster in the second quadrant (see below), 
(iii) Although the polar motion excitations by the latest 6 events have magnitude 
that were hardly detectable. Events I and II, if happened today, could be readily 
detected in the polar motion measurement. It should be noted, however, that it is 
in the polar motion excitation that an earthquake leaves its signature as a step 
function. In the polar motion path itself (which represents a time integration of the 
excitation, see equation ( 1)), the earthquake manifests itself as a “kink” in the path, 
and hence more difficult to detect. Further difficulties arise because of the abrupt 
nature of tbe earthquake source and the relatively “gappy” sampling of earth 
rotation in the past. Only recently has a virtually continuous monitoring of earth 
rotation become feasible by means of the GPS technique. 

Figure 3 compares the power spectra of the observed and the seismic excitations 
for the same period of 1977.0-1993.6. The spectra are computed using Thomson s 
( 1982) multitaper technique after removal of the mean values. The different spectral 
characteristic is evident. The frequency dependence/" of the spectrum of the seismic- 
excitation is found to follow n = 2.0, similar to a Brownian motion process. The 
power of the seismic T is in general 40 60dB lower than that observed. In 
particular, the power difference at the Chandler frequency band around 0.83 cycle 
per year is about 45 dB. We have also calculated and examined possible correlation 
between the two excitation functions in the Chandler band. Only a weak coherence 
at a moderate confidence level was found. These findings are consistent with Gross 
(1986). 

Despite the small magnitude, the seismic excitations T collectively are interest- 
ing in their own right. Following Chao and Gross (1987), we shall now conduct 
statistical analyses to reveal their peculiar behavior. 

Figure 4 shows the angular histogram (“rose diagram”) of the arguments of the 
1 1,015 T’s in thirty-six 10 increment bins. Apart from a concentration around 
15 E, an abnormally large number of 'P’s cluster around 140 'E. The distribution 
pattern is remarkably similar to the 2,146 T’s analyzed by Chao and Gross 
(1987), indicating that this pattern is robust with respect to time and number of 
earthquake samples. In fact, the statistical tendency of this angular anomaly is 
stronger now with the substantially greater number of samples: the normalized x’ 
found here is as high as 14.1, compared to 3.81 found in Chao and Gross ( 1987). 
Compared with, say, the V'A> significant level of 1.64 or the 0.1% significant level of 
1 .90 for a random distribution (at 35 degrees of freedom), this asserts the extremely 
non-random nature of the distribution of Figure 4. 
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Figure 3 

The ( multitapered) power spectrum (in dB units of mas^/cpy) of the observed polar-motion excitation 
compared with that induced by earthquakes during 1977.0 1993.6. Mean values have been removed 

from the time series. 


The preference of earthquakes in nudging the rotation pole towards ~ 140"E is 
also evident in Figure 2. The straight lines are the least-squares fit to the curves; 
their slopes give the velocity of the polar drift induced by earthquakes. This polar 
drift velocity vector (0.019 mas/yr, 131 "E), as plotted in Figure 5, indeed points to 
that general direction. The magnitude of the vector is here magnified by 100 times, 
in order to be shown against other vectors: 

The one labeled '‘O&D’' (4.5mas/yr, 148 E) is that similarly obtained using 
O’Connell and DziewonskTs (1976) calculation for 30 great earthquakes during 
1900-1964. So it appears that the earthquake-induced polar drift has continued its 
journey toward ^ 140°E at least since 1900 when fairly reliable seismic records had 
become available. It should be cautioned, however, that O’Connell and 
Dziewonski may very well have greatly, and presumably to different extent, 
overestimated the sizes of their studied earthquakes (Kanamori, 1976). Even 
though the directions of the polar excitation by individual earthquakes, which are 
determined by the source mechanisms according to plate tectonics, are realistic, the 
(vectorial) cumulative drift direction may be biased. However, closer examination 
of these individual directions indicates that this bias is not severe: Figure 6 plots the 
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Figure 4 

The angular histogram (^Tosc diagram") for the direction of the seismic excitation of polar motion by 
major earthquakes for 1977.0 1993.6, in 36 angular bins with respect to the terrestrial coordinate 

system. 


rose diagram (in 20 bins) of O'Connell and Dziewonski's 30 polar excitation 
directions and their cumulative polar drift. The directions have the same heavy 
clustering in the second quadrant, and hence the cumulative drift is largely 
unidirectional. 

The other two vectors in Figure 5 are from polar motion observations; they 
agree well with each other. The one labeled “Space93” (3.84mas/yr, — 68 E) is 
simply obtained from the fitted line in Figure 1. The one labeled “Pole93” 
(3.22mas/yr, -81 E) is similarly obtained from a polar motion dataset for 
1900-1993, primarily based on the International Latitude Service data since 1900. 
It is seen that the observed polar drift directions are roughly opposite to those 
induced by earthquakes. 


4. Conclusions 

The direction of the seismic excitation depends on the focal location and source 
mechanism. Through computation, it is found that, at least since 1900, the seismic 
excitations collectively have a statistically strong tendency to nudge the pole 
towards ~ 140 E, away from the actually observed polar drift direction. This non- 
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Figure 5 

The calculated polar drift velocity induced by major earthquakes for 1977.0-1993.6, magnified 100 
times, in comparison with O'Connell and Dziewonski's estimate (labeled “O&D") and the observed 
polar drift (labeled ‘"Spacc93" and “Pole93"). 


random behavior is similar to other earthquake-induced changes computed for 
length-of-day and low-degree gravitational field, which indicate a strong tendency 
for the earthquakes to make the earth rounder and more compact (Chao and 
Gross, 1987). Specific questions arise: Why do earthquakes seem to ‘‘recognize'’ 
the existence of the rotation axis? Is there any as yet unseen, dynamic connection 
between the pole position and the occurrence and source mechanism of earth- 
quakes? Or are they manifestations of some behind-the-scene geophyiscal processes? 
The dynamic reasons for such peculiar behavior of earthquakes are not clear in 
terms of the grand scheme of the plate tectonics. At the present time these remain 
mere speculations. 

The magnitude of individual earthquake effect, on the other hand, depends 
largely on the seismic moment of the event. The rule of thumb is that the seismic 
moment of 10“^ Nm would roughly produce 1 mas in polar-motion excitation and 
1 /ts in length-of-day change; the actual values depend on the focal location and 
seismic source mechanism. The polar motion excitations produced by the largest 
earthquakes during 1977-1993 were still an order of magnitude too small to be 
detected even with today's measuring accuracy of about 1 mas. However, it was 
calculated that the 1960 Chile event should have produced a discontinuity as large 
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Figure 6 

(a) The rose diagram for the direction of the seismic excitation of polar motion by 30 great earthquakes 
during 1900 1964 studied by O’Connell and Dziewonski, in 18 angular bins, (b) Their cumulative polar 
excitation in the terrestrial coordinate system. 

as 23 mas in the polar-motion excitation function (but only 8 //s in ALOD). So an 
earthquake of only one tenth the size of that event, if happened today, could be 
comfortably detected in polar-motion observations, especially now that the GPS 
technique routinely provides sub-daily temporal resolutions. Further improvements 
in observation should allow detection of smaller (and hence potentially more 
numerous) earthquake signatures in the future. 
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